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Abstract 

In this paper, we categorize non-rotating black hole spacetimes based on their pole 
structure and in each of these categories we determine whether the WKB approxi- 
mation is a valid approximation for calculating the asymptotic quasinormal modes. 
We show that Schwarzschild black holes with the Gauss-Bonnet correction belong to 
the category in which the WKB approximation is invalid for calculating these modes. 
In this context, we further discuss and clarify some of the ambiguity in the literature 
surrounding the validity conditions provided for the WKB approximation. 



1. Introduction 



The equations governing various classes of non-rotating black hole quasinormal mode 
(QNM) perturbations in radial coordinate r have the general form 

d 2 ^ 

— + R(r)* = 0, (1.1) 
with the function R(r) of the general form 



^( r ) = JTT2 ^ - V ^ 
f{r) 2 1 
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where u is the complex QNM frequency with real (ojr) and imaginary (ojj) parts, V(r) is 
a function that depends on the type of perturbation and the spacetime structure, and f(r) 
is the black hole metric function that appears in the spcetime metric 

dv^ 1 

ds * = -/(r)rft 2 + _ + r 2 dnl_ 2 (1.3) 
f{r) 

with the spacetime dimension D. 

Methods developed to calculate the highly damped [U El |3] and other asymptotic [U [5] 
QNM frequencies of black holes (see also [6] for a comprehensive review of black hole 
QNM's) are based on the WKB approximation. Here, by highly damped we mean 
\u\ ~ \u)i\ » ujr and by asymptotic we mean when \u\ — > oo in general. The highly 
damped QNM's have attracted considerable attention in the literature as a result of the 
work of Hod [7] and more recently Maggiore|8j, who have attempted to link these modes 
to the horizon area quantization of black holes. Also, in a recent paper [9], a link be- 
tween the highly damped QNM's and the small scale structure of black hole spacetimes 
has been established. The authors of [9] show that the highly damped QNM's are sensitive 
to any small scale (for example a small charge or angular momentum) that is added to 
a Schwarzschild metric. This includes any additional length scale due to a quantum cor- 
rection. The authors of [9] demonstrate this explicitly by calculating the highly damped 
QNM's of the quantum corrected black hole derived in [10]. The results of this paper are 
also relevant to the analytic calculations of the greybody factors for black holes at large 
imaginary frequencies. See, for example, |llj . 

The two approximate, linearly independent, complex, WKB solutions to the differential 
equation (11. ip are 

#W( r ) = Q-V2( r ) exp[+2 f t Q{r')dr'} 

(1.4) 

#W( r ) = Q-V2( r ) exp [-z f t r Q(r')dr'\ , 

where Q 2 = R and t is customarily taken to be a zero of Q in the complex r-plane. In 
the methods developed in [H El 131 HI E], one needs to analytically continue the complex 
WKB solutions (ll.4jl into the whole complex coordinate plane. After locating the poles and 
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zeros of the function R(r), and consequently Q(r), one is able to determine the topology of 
certain lines called Stokes and anti-Stokes lines. Along Stokes lines the phase of the WKB 
solutions (/ Qdr) is purely imaginary and along anti-Stokes lines this phase is purely real. 
In order to determine a WKB condition on the QNM frequency u, in most cases one needs 
to find a closed contour around the pole at the event horizon. We choose this contour along 
anti-Stokes lines because the solutions do not change character on these lines. In other 
words, if we know the behavior of the solution on a point of an anti-Stokes line, then we 
know the behavior of the solution everywhere on that line. Along this contour, one always 
encounters some of the poles of the function R(r) in the complex planeQ The details of 
these calculations are not important for the purpose of this paper. What is important is 
that, in the infinite damping limit, it is necessary to approach infinitely close to some of 
these poles. This can be realized by noticing that in the limit \u>\ — > oo, the function V(r) 
can be neglected everywhere in the complex plane compared to u> 2 except when we are 
infinitely close to the poles of V(r). Consequently, one needs to make sure that the WKB 
approximation is valid in the region near those poles. 

The rest of this paper proceeds as follows. In Sec. 2, we provide a generic function 
R(r) which includes the pole structure of all types of non-rotating black hole spacetimes 
(at least the ones we are aware of). In Sec. 3, we categorize black hole spacetimes into four 
categories and in each of these categories we determine whether the WKB approximation 
is valid for calculating the asymptotic QNM frequencies. In Sec. 4, we show that the WKB 
approximation is not valid for calculating the highly damped QNM's of Schwarzschild black 
holes with the Gauss-Bonnet correction. In Sec. 5, we discuss additional WKB validity 
conditions in the context of the black hole QNM problem. 



2. A Generic Example 



Let us assume that R has the generic form 

*(0 = ^ - j^ry a . (2-i) 

where r* is the location of one of the poles of the function R(r) in the complex plane, A is 
a real or complex constant, and a is a real number greater than zero. We assume that the 
function F(r) is a general function that can have a pole at the same location of r — r* but 
with a power less than a in its denominator. We want to check if the WKB approximation 
is valid in the region where r — > r*. A well known condition for checking the validity of 
the WKB approximation is: 

, „ 3 dR . 

iT*-H « 1 • (2.2) 
dr 

There exist some ambiguity in the literature regarding the extent of the applicability of 
the above condition. Many books [121 1131 HH E3 interpret ( 12. 2 p as a condition that the 



In the method developed in [T], one does not need to determine the topology of Stokes lines or the 
position of the zeros of Q. For more details see the end of Sec. 4. 
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function R be "slowly varying" and suggest it is sufficient to guarantee accuracy of the WKB 
approximation. However, in (16] this condition is only one of two conditions required for 
sufficiency. Perhaps more accurately, condition (I2.2p is referred to as a "rule of thumb" 
in [IT]. In this section, we will discuss the extent of the validity of this condition in the 
context of the generic example provided in Eq. (12. f p . We will discuss other WKB validity 
conditions in Sec. 5. 

For a > 2, it is easy to show that near r = r* the validity condition in (I2.2p is satisfied. 
In order to further make sure that the WKB approximation is valid for a > 2, we need to 
check if, in the asymptotic region r — > r*, the solution to the wave equation (11. ip matches 
the WKB solutions (11.41) . Therefore, we solve the wave equation (11.11) with the function 
R(r) given in (12.11) in the limit r — > r*. Note that, in this limit, F(r) can be neglected. 
The solution is 



2-a 



V— A(r — r*) 



2-a 
2 



+C 2 \/r — r* J__j_ 



2-a 



V—A(r - r*) 2 



(2-3) 

where are constants and J v is the Bessel function of the first kind. On the other hand, 
the WKB solutions (II .4p behave like 



^wkb{t) « C+(r -r*)Va- 



-A(r-r. t )~ 



+ C-(r 



rJ*e 



-A(r— r»)~ 



(2.4) 



when r r*. For the WKB approximation to be valid near r = r*, the solutions (I2.3P and 
(12 .4p have to match. To check this, we use the asymptotic behavior of the Bessel function 
J,Ax) when x » 1, i.e. 




J,Ax) «( — 



2 / w i 

cos x 

ttx V 2 4 



(2.5) 



Using the above asymptotic behavior, we can easily show that solution (12. 3p can be reduced 
to the WKB solution (12. 4p . Here, we have to point out that the asymptotic behavior (12. 5 p 
can only be used when the argument of the Bessel function has a real value. Note that 
the argument of the Bessel function in (12. 3 p is equal to the phase (/ Qdr) of the WKB 
solutions. Since at every point in the complex r-plane we can find a curve on which / Qdr 
is purely real, we can always use the asymptotic behavior (12. 5p . 



Things change for the case of a 



2. In this case, when r r„, we have 







2 


R~l dR 






dr 







(2.6) 



which tells us that the condition (I2.2p is only satisfied if \A\ » 1. However, we can still 
use the WKB approximation if we apply a modification. It is easy to show that, when 
rRir*, the solutions of the wave equation (II. ip with R given by (12.11) are 



^f(r) oc (r — r* 



V1+4A 
l 2 2 



Meanwhile, if we take Q 2 = R, we get 



Q 2 oc ( r - r*)2 



(2.7) 



(2- 
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and 

J Qdr « ±iVJ\n(r - r*) , (2.9) 

which means the WKB solutions are 

*i,a oc (r-r*)^ ±V3 . (2.10) 

Therefore, the WKB solutions do not have the correct behavior near r = r*. This discrep- 
ancy, however, can be resolved by choosing 

Q 2 = i?~ A , 1 , 2 ■ (2.H) 
4(r — r*)^ 

This modification is very similar to the Langer modification used in the WKB approxima- 
tion of radial quantum problems. It is important to point out that the above shift does 
not change the topology of anti-Stokes lines, which are crucial in the asymptotic QNM 
frequency calculations. 

In the case of a < 2, condition (I2.2p is not satisfied near r = r*. The fact that the 
WKB approximation is not valid for a < 2 can be observed by noticing that the solution 
( 12. 3 p approaches 

tt(r) ^ C x {r - r*) + C 2 (2.12) 
in the r — )■ r* limit. But, in the same limit, the WKB solution ( 12. 4 p approaches 

*WKB(r) « C + (r - r,)t + C_(r - r*)t . (2.13) 

It is clear that the WKB solutions do not have the correct behavior near these poles and 
their behavior cannot be corrected by a simple shift of the form (12. lip in which the topology 
of anti-Stokes lines is unaltered. To be more precise, in order to match the solutions of the 
QNM wave equation and the WKB solutions near r = r*, we need to take 

Q 2 = ~, ~~ T5" (2.14) 

(r — r*)2 

if the constant Ci is zero. Otherwise, we need the shift 

Q 2 = R{r-r^ . (2.15) 

Such shifts will modify the topology of Stokes/anti-Stokes lines since the pole structure is 
altered. 

This example shows that if one accounts for the Langer-type modification given in 
(12. lip , then the WKB validity condition ( 12. 2 j) is a sufficient condition. However, if this 
validity condition is not satisfied, it may still be possible to use the WKB approximation. 
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3. Categorizing Black Hole Spacetimes for QNM Calculations 



Based on the generic function R(r) provided in Eq. (12. ip . we can categorize non-rotating 
black holes based on their pole structure into four different categories: 

1) At all poles r*, a > 2. 

2) At all poles r*, a = 2. 

3) At all poles r*, a < 2. 

4) Some of the poles have a > 2 and others have a < 2. 

Some black hole examples of the second category are Schwarzschild, Reissner-Nordstrom, 
Schwarzschild-de Sitter, Schwarzschild-anti de Sitter and the non-singular quantum cor- 
rected black hole spacetime derived in [10] . We are not aware of any example of the first or 
the third category. The Schwarzschild black hole with the Gauss-Bonnet correction, which 
will be explained in the next section, is an example of the fourth category. The WKB 
approximation in calculating the asymptotic QNM frequency is always valid for the first 
and second category and it is always invalid for the third category . The validity of the 
WKB approximation in the fourth category depends on whether we encounter a pole with 
a < 2 along the anti-Stokes lines encircling the event horizon or not. If we do, then the 
WKB approximation is invalid. Otherwise, it is still fine to use this approximation. 

In this paper we have only talked about non-rotating black holes. For rotating black 
holes, similar categorization may be possible. For example, the pole structure of Kerr black 
holes is of category 2, but in the case of Kerr black holes the asymptotic behavior of the 
wave equation near the poles is irrelevant in the high damping limit because the contour 
that we take along anti-Stokes lines never encounters any pole in the complex plane [T%] . 
Therefore, in the Kerr case, we do not need to check if the WKB approximation is valid 
near the poles. We predict that the same should be true for other types of rotating black 
holes such as Kerr-de Sitter, Kerr-anti de Sitter, etc. 



4. The Case of Gauss-Bonnet 



Schwarzschild black holes with Gauss-Bonnet (G-B) corrections [T^J [213 EI] (hereafter 
G-B black holes) are black holes in spacetime dimensions greater than four which appear 
when one includes the leading order higher curvature terms arising in the low energy limit 
of string theories [22]. 

The authors of [23] use a combination of analytic and numeric techniques based on 
the WKB approximation to calculate the QNM frequencies of G-B black holes in the 
intermediate damping range of 



1 << \u\jJL 1 



« 



4(L>-2) 2 -l 



(8a/i) D - 



2D-_4 
1 



D- 



(4.1) 



where D is the spacetime dimension, a is the G-B coupling constant and the parameter \i 
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is related to the ADM mass M given by 



M = . (4.2) 



Working in the intermediate damping region mentioned above allows the authors of 

to do the calculations in a region of the complex r-plane where r >> (Safi) - 1 . In 
this region, the calculations are not affected by the fine structure of the extra poles in the 
complex r-plane, which appear as a result of the G-B correction. These poles are located at 
the roots of r D_1 + 8a/x = and (D — 3)r D_1 +Aafi(D — 5) = 0. Under such circumstances 
it was shown in [23] that, in the intermediate damping region (14. ip . 

2tt 2 ^ + a cor « ln(2.99926) (4.3) 
y2yU — a 

when D = 5. This result is very close to the Schwarzschild result of ln(3)[lj. In the limit 
of a — > 0, the correct Schwarzschild value of ln(3) is recovered. For further discussions on 
the Schwarzschild limit of the highly damped QNM's of different black holes, see [2] and 
[SJ. In the highly damped limit, where 



M >> 



A(D - 2) 2 - 1 

/ , 2D-A 

(8a/x) v- 1 



// , (4.4) 



the authors of [23] found that the standard analytic techniques cannot be applied in a 
straightforward way. In what follows we show that, in fact, the WKB approximation is an 
invalid approximation for calculating the QNM's of G-B black holes in the highly damped 
region (14.4H . 

We can use the metric function for a Schwarzschild black hole with the G-B 
correction |19] 

/W = 1+ 2^-^V 1 + ^ (45) 

to find that the function R(r) in five spacetime dimensions can be approximated in the 
\u)\ —> oo limit by [23] 

J 2 35 48a/i 

R y r > + ( r 4 + 8a/x )3/ 2 • t 4 - b J 

This function has a pole structure of category 4, where for one pole a = 2 and for the 
other poles a = 3/2. Therefore, according to the results obtained in Sec. 2, the WKB 
approximation would be invalid if one encounters a pole with a = 3/2 along anti-Stokes 
lines that loop around the event horizon. This turns out to be the case for the G-B black 
hole QNM's in the highly damped region (I4.4p . The path along anti-Stokes lines for the G- 
B case in the high damping limit is shown in Figure 1 of reference [23J . In the intermediate 
damping region (14. ip . however, the path that we take along anti-Stokes lines around the 
event horizon never passes near the poles at r = (— 8a/i) 1//4 . Therefore, the calculations 
can be done in the intermediate damping region using the WKB approximation. The path 
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along anti-Stokes lines for the G-B case in the intermediate damping region can be found 
in Figure 4 of reference [23J . 

We also would like to point out that in the method developed by Motl and NeitzkepQ 
for calculating the highly damped QNM's, determining the location of the zeros of Q or 
the fine structure of Stokes and anti-Stokes lines near the poles is not necessary. In this 
method, the wave equation near the poles is solved analytically and then the solutions 
are matched to the WKB solutions along anti-Stokes lines when one moves away from 
the poles. In the region close to the poles where r 4 — > — 8a/i the wave equation can be 
approximated as 



d 2 ^ 

dr 2 



~2/7 



48a/i 



+ 8a/x) 3 / 2 



^ = 



(4.7) 



The WKB validity condition in ( 12. 2 p is not satisfied for the above wave equation because 



dr 



2/x 



+ 



48a/i 



-3/2 



-288a// 



(r 4 + 8a/x) 5 /2 



oc 



y + Sa/i)" 1 / 4 > 1 (4 



near the poles at r = (8a/i) 1//4 . Unfortunately, Eq. (14. 7p does not have a simple solution 
and we cannot show directly if the solutions to this differential equation can be matched 
to the WKB solutions ( II. 4p . but the fact that the validity condition ( 12. 2 j) is not satisfied 
is a strong indication that the WKB approximation is invalid. In addition, since the two 
methods in pQ and [2] should produce consistent results, it would be highly unlikely that 
in one method the WKB approximation is valid and in the other invalid. 



5. Other Validity Conditions for the WKB Approximation 



To discuss some of the other validity conditions used for the WKB approximation in 
the literature, it is useful to review the derivation below. For more details see [24J. 

In the WKB theory, we assume a solution to the wave equation (II .ip in terms of an 
exponential power series of the form 



\I/(r) ~ exp 



1 oo 
6 n=0 



(5.1) 



when 5—t-O. The WKB series Y^3 n S n { r ) is divergent for most cases and for that reason 
the asymptotic notation ~ is used in ( 15. ip instead of the equality sign (=). By inserting 
(15. ip into the wave equation (11.11) . we obtain 



±S» + - 5 S' S[ + ±£% + ... = -R(r) 



(5.2) 



The first term on the left is the dominant term as 5 — > and must have the same order of 
magnitude as R{r). Therefore, we get 

1 



Q 



r2 



-R(r) 



(5.3) 
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The coefficients of 5 l , i > —2, in ( 15. 2 p have to vanish. This leads to 

2S' S[ + 5* = , 

and 

n-l 

2S' S' n + S'U + = , n > 2 . 

3=1 

Solving these equations gives 



5 



So(r) 



± 



(5.4) 
(5.5) 

(5.6) 
(5.7) 

-Rfl 2 ' 32(-i?) 5 /2_ ' (5 - 8) 

and so on. The WKB approximation used in ( II. 4p is just exp[|So + Si], which is sometimes 
called the leading-order WKB approximation. 

For the WKB approximation to be valid on an interval, the series ^Z5 n S n {r) has to be 
an asymptotic series in 5 as 5 — > uniformly for all r on the interval. This requires that 

:S Q (r) as 5^0 



Si(r) 



5So(r) 



± 



1 

~4 
R" 



-R(t)dt 
H-R(r)} , 



+ 



5R 



r> 



dt 



StM « 



5 



(5.9) 



and 



\5S 2 (r)\ « |^i(r)] as 5^0, (5.10) 

which have to hold uniformly in r. We also need to require that the first term we drop is 
small, i.e. 

\SS 2 (r)\ « 1 as 5^0. (5.11) 

This is because the WKB series appear in the exponent and if the above condition is not 
satisfied the error can be large. Here we are only concerned with the leading-order WKB 
approximation, which means that the conditions can stop at S 2 and no further conditions 
are necessary. 

The validity condition (12. 2p . which was mentioned in Sec. 2, implies the condition 
( 15. lip . This can be shown by rewriting Eq. ( 15. 8 p as 



SS 2 (r) = T 



R' 



8(-R) 3 / 2 



± 



R 



12 



32(-i?) 5 /2 



dt . 



(5.12) 



It is now easy to see that \5S 2 \ « 1 as long as condition (12.21) is satisfied. It is also 
possible to see where (12.21) comes from by realizing that the dominant term, Sq/S 2 , in Eq. 
(15. 2p has to be much larger than either one of the terms of O as 5 — > 0. 
It is interesting to point out that, similar to the validity condition 



the sets of 

validity conditions ( 15. 9p through ( 15. lip are all satisfied for category 1 black holes and fail 
for category 3. For category 2 black holes, they may or may not be satisfied depending 
on how large or small the constant A is in Eq. (II. 2p . In such a situation, a Langer-type 
modification will allow us to use the WKB approximation. 
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